A one dimensional analysis of turbulence and 
its intermittence for the d-dimensional 
stochastic Burgers equation 

o 

O . AD Neate A Truman 

(N ■ 

I Department of Mathematics, University of Wales Swansea, 

3 ■ Singleton Park, Swansea, SA2 8PP, Wales, UK. 



00 



ON 

in 



X 

5-H 



February 1, 2008 



Abstract 

i -pH ■ The inviscid limit of the stochastic Burgers equation is discussed 

in terms of the level surfaces of the minimising Hamilton-Jacobi func- 
tion, the classical mechanical caustic and the Maxwell set and their 
algebraic pre-images under the classical mechanical flow map. The 
problem is analysed in terms of a reduced (one dimensional) action 
^ \ function. We demonstrate that the geometry of the caustic, level sur- 

faces and Maxwell set can change infinitely rapidly causing turbulent 
behaviour which is stochastic in nature. The intermittence of this tur- 
bulence is demonstrated in terms of the recurrence of two processes. 



o 

O '. 1 Introduction 



Burgers equation has been used in studying turbulence and in modelling the 
large scale structure of the universe [U EH EE], as well as to obtain detailed 
asymptotics for stochastic Schrodinger and heat equations. pUi EH EE1 EH 
EDI El]. It has also played a part in Arnol'd's work on caustics and Maslov's 
works in semiclassical quantum mechanics [31 IU EE ED] ■ 

We consider the stochastic viscous Burgers equation for the velocity field 
v>*(x, t) E R d , where x E R d and t > 0, 

— + («" • V) if = y Av " ~ W(x) - eVk t {x)W t , 

with initial condition 0) = VS (x) + 0(/i 2 ). Here W t denotes white 
noise and /i 2 is the coefficient of viscosity which we assume to be small. 
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We are interested in the advent of discontinuities in the inviscid limit 
of the Burgers fluid velocity v°(x,t) where v ,J, (x, t) -> v°(x,t) as fi —> 0. 
Using the Hopf-Cole transformation v^(x,t) = -/i 2 Vlnti''(i,t), the Burgers 
equation becomes the Stratonovich heat equation, 

- ^-Au^ + ^- 2 V(x)u" + -zktWu 1 * o W t , 



dt 2 v ' fi 

with initial condition u^(x, 0) = exp f~ ~^ l ^o0e)> where the convergence 
factor To is related to the initial Burgers fluid density. 
Now let, 



A{X}:=\f X 2 (s)ds- [ V(X(s))ds-e[ k s (X(s))dW s , (1) 
^ Jo Jo Jo 

and select a path X which minimises This requires, 

dX(s) + VV(X(s)) ds + eVk s (X(s)) dW s = 0. 

We then define the stochastic action A(X(0),x,t) := inf {^4[X] : X(t) = x} . 
Setting, 

A(X(0),x, t) := S (X(0)) + A{X{0), x, t), 

and then minimising A over X(0), gives X(0) = VSq(X(0)). Moreover, it 
follows that, 

S t (x) := inf {A{X(0),x,t)}, 

X(0) 

is the minimal solution of the Hamilton- Jacobi equation, 

dS t + Q| V^l 2 + V(x)^j dt + ek t (x) dW t = 0, S t=0 (x) = S (x). (2) 

Following the work of Donsker, and Freidlin and Wentzell [12] , — /i 2 In ^(x, t) - 
St(x) as \i — > 0. This gives the inviscid limit of the minimal entropy solution 
of Burgers equation as v°(x, t) = VS t (x) [5]. 

Definition 1.1. The stochastic wavefront at time t is defined to be the set, 

W t = {x: S t (x) = 0}. 

For small /i and fixed t, t) switches continuously from being expo- 
nentially large to small as x crosses the wavefront Wt- However, and w M 
can also switch discontinuously. 
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Define the classical flow map $ s : M. d — > M. d by, 

d$ s + W($ s ) ds + eVk s ($ s ) dW s = 0, $ = id, $ = V5 . 

Since X(t) = x it follows that A(s) = $ s where the pre-image 

a; (a;, t) = $^ (x) is not necessarily unique. 

Given some regularity and boundedness, the global inverse function the- 
orem gives a caustic time T{uj) such that for < t < T(u), $ t is a random 
diffeomorphism; before the caustic time v°(x, t) = $ t ($ < ~ 1 (x)) is the inviscid 
limit of a classical solution of the Burgers equation with probability one. 

The method of characteristics suggests that discontinuities in v°(x, t) are 
associated with the non-uniqueness of the real pre-image xq(x, t). When this 
occurs, the classical flow map $i focusses an infinitesimal volume of points 
dxo into a zero volume dX(t). 

Definition 1.2. The caustic at time t is defined to be the set, 



Assume that x has n real pre-images, 

^t 1 i x ) = {%o(l)(x,t),x (2)(x,t), . . . ,x (n)(x,t)} , 

where each xo(i)(x,t) e M. d . Then the Feynman-Kac formula and Laplace's 
method in infinite dimensions give for a non-degenerate critical point, 



,t)=y* J 9 i exp ( 
i=i \ 



u^t) = yj iew [-^^), (3) 
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where Sq(x, t) := S (x (i)(x, t))+A (x (i)(x, t),x, t) , and 0, is an asymptotic 
series in /i 2 . An asymptotic series in /i 2 can also be found for v^(x,t) [32] . 
Note that S t (x) = minjS'^x, t) : i — 1, 2, . . . , n}. 

Definition 1.3. TTie Hamilton- Jacobi level surface is the set, 

HI = {x : Sq(x, t) = c for some i\ . 
The zero level surface H® includes the wavefront W t . 

As \x — > 0, the dominant term in the expansion comes from the min- 
imising x (i)(x, t) which we denote x (x,t). Assuming x (x, t) is unique, we 
obtain the inviscid limit of the Burgers fluid velocity v °(x, i) = <& t (xq(x, t)) . 
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If the minimising pre-image Xq(x, t) suddenly changes value between two 
pre-images x (i)(x,t) and x (j)(x,t), a jump discontinuity will occur in 
v°(x, t), the inviscid limit of the Burgers fluid velocity. There are two distinct 
ways in which the minimiser can change; either two pre-images coalesce and 
disappear (become complex), or the minimiser switches between two pre- 
images at the same action value. The first of these occurs as x crosses the 
caustic and when the minimiser disappears the caustic is said to be cool. The 
second occurs as x crosses the Maxwell set and again, when the minimiser is 
involved the Maxwell set is said to be cool. 

Definition 1.4. The Maxwell set is, 

M t = {x : Bx ,x e M. d s.t. 

x = $t(x ) = ®t{xo), x ^ x Q and A(x , x, t) = A(x Q , x, t)} . 

Example 1.5 (The generic Cusp). Let V(x, y) = 0, k t (x, y) = and So(xo, 2/°) = 
Xgt/o/2. This initial condition leads to the generic Cusp, a semicubical 
parabolic caustic shown in Figure [H The caustic C t (long dash) is given 





and the Maxwell set M t (short dash) is x = for y > —1/t. 




Figure 1: Cusp and Tricorn. 
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Notation: Throughout this paper x, x , x t: xu c ) etc will denote vectors, where 
normally x = $> t (xo). Cartesian coordinates of these will be indicated 
using a sub /superscript where relevant; thus x = {x\, X2, ■ ■ ■ , xj), xo = 
(xj, Xq, . . . , xf) etc. The only exception will be in discussions of explicit 
examples in two and three dimensions when we will use (x,y) and (xo,yo) 
etc to denote the vectors. 



2 Some background 

We begin by summarising some results of Davies, Truman and Zhao (DTZ) 
[HlEj. Following equation ([T]), let the stochastic action be defined, 



A{x ,p ,t) = \ f X{s) 2 ds- [ 
1 Jo Jo 



t r 

V(X(s)) ds + ek s (X(s)) dW s 



where X(s) = X(s,x ,p ) G M. d and, 

dX(s) = -VV(X(s)) ds - eVk s (X(s)) dW s , X(0) = x , X(0) = p , 

for s G [0,t] with x ,p G M. d . We assume X(s) is T s measurable and unique 
where T s is the sigma algebra generated by X{u) up to time s. It follows 
from Kunita 1181: 



Lemma 2.1. Assume Sq,V G C 2 and k t G C 2,0 , W,Vk t Lipschitz with 
Hessians V 2 V, V 2 k t and all second derivatives with respect to space variables 
of V and k t bounded. Then for p , possibly xq dependent, 

^(x ,po,t) = X(t)-?^-X a (0), a=l,2,...,d. 

Methods of Kolokoltsov et al [17] guarantee that for small t the map 
p i— > X(t,x ,po) is onto for all x . 

Therefore, we can define A(x ,x,t) := A(x ,p (x , x,t),t) where p = 
Po(xo,x,t) is the random minimiser (which we assume to be unique) of 
A(xo,po,t) when X(t,xo,po) = x. Thus, the stochastic action corresponding 
to the initial momentum VSo(x ) is <A(xo, x,t) := A(x ,x, t) + S (x ). 

Theorem 2.2. If $t is the stochastic flow map then, 

d 

$t(x ) = x 7— [A(x ,x,t)} = 0, a = l,2,...,d. 

ox 
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The Hamilton- Jacobi level surface H£ is found by eliminating xo between, 

dA 

A(x ,x,t) = c, — (x ,x,t) = a = l,2,...,d. 

Alternatively, if we eliminate x to give an expression in xq, we have the pre- 
level surface ^ 1 H^. Similarly the caustic Ct (and pre-caustic $^~ 1 C't) are 
obtained by eliminating xq (or x) between, 

/ Q2j{ \ dA 

det [-^J^ x ^) =0 ' ^' I ' i ) = a = l ) 2,...,d. 

VCXqCXq ' a ,p=i,2,...,d 

These pre-images are calculated algebraically and are not necessarily the 
topological inverse images of the surfaces C t and if t c under 

The caustic surface can be parameterised using its pre-image by ap- 
plying the stochastic flow map (a pre-parameterisation) . This allows us 
to control the domain of the pre-images and in particular restrict them 
to real values. If we can locally solve the equation of the pre-caustic to 
give, xl = Ai, Xq = A 2 , . . . , x 1 = Xd-i, Xq = Xq c (Ai, A 2 , . . . , Xd-i) , 
then the pre-parameterisation of the caustic is x t (A) := $ t (A,Xq C (A)) where 
A = (A 1 ,...,A d _ 1 )eM d - 1 

We next outline a one dimensional analysis first described by Reynolds, 
Truman and Williams (RTW) 



Definition 2.3. The d- dimensional flow map $ t is globally reducible if for 
any x = (xi, x 2 , . . . , xj) and x = (xq, Xq, . . . , Xq) where x = $ t (x ), it is 
possible to write each coordinate Xq as a function of the lower coordinates. 
That is, 

x = $t(xo) x o= Xq(x, Xq, Xq, . . . , Xq -1 , t) for a — d, d — 1, . . . , 2. 

(4) 

Therefore, using Theorem 12.21 the flow map is globally reducible if we 
can find a chain of C 2 functions Xq, Xq -1 , . . . , Xq such that, 

X = X 0\ X i ^0' ^0' ■ ■ - x o yt) -vv- ^ (xo, x, t) = 0, 

-1 -1 12 d 2 &4 ] 2 d 

"^0 = "^0 ( X ' X Qi X Qi ■ ■ ■ x yt) -vv- - (Xq, Xq, . . . , Xq(. . .), X, t) — 0, 

<7Xq 
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dA 



(Xq, Xg, Xq(x, Xq, Xq, t), . . . , Xq(. . .), X, t) — 0. 



This requires that no roots are repeated to ensure that none of the second 
derivatives of A vanish. We assume also that there is a favoured ordering of 
coordinates and a corresponding decomposition of $ f which allows the non- 
uniqueness to be reduced to the level of the x\ coordinate. This assumption 
appears to be quite restrictive. However, local reducibility at x follows from 
the implicit function theorem and some mild assumptions on the derivatives 
oiA (see [22]). 



Definition 2.4. IfQt is globally reducible then the reduced action function is 
the univariate function obtained by evaluating the action with equations (fjjj, 

f(x,t){ x o) := f{ x 0i x i ^) = ^(-^O' x o( x i x Qi ^)) x o(- ■■)■>■■■, x,t). 
Lemma 2.5. If &t is globally reducible, modulo the above assumptions, 



det 



3 2 A 

(dx ) 2 



(x ,x,t) 



n 

a=l 



I Air 1 r a T a+1 ( 

8xq) °' 1 



.),x,t) 



r 0~ x o( x i x O^) 



where the first term is f" xt ){xl) and the last d — 1 terms are nonzero. 

Theorem 2.6. Let the classical mechanical flow map $t be globally reducible. 
Then: 

1. f[ xt ^{x\) = and the equations (f^Jj <^> x = $ f (xo), 

2 - f(x,t)( x o) = f(x,t)( x o) = and the equations (fj]j 

x = $t(^o) is such that the number of real solutions xq changes. 

Corollary 2.7. Let x t (X) denote the pre-parameterisation of the caustic with 
A = (Ai, A 2 , • • • , X d -x) e R*" 1 . Then f ixtm) {Xi) = f^ m) (^) = 0- 

Consider an example where for x on one side of the caustic there are 
four real critical points on f( Xt t){x\) enumerated xl(i)(x, t) for i — 1 to 4, 
and denote the minimising critical point by xl(x,t). Figure [2] illustrates how 
the minimiser jumps from (a) to (b) as x crosses the caustic if the point of 
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Before Caustic 

xj(3) 




xl(2)=xl(x,t) 



Minimiser at 
x\{2){x,t) = xl{x,t) 



On Cool Caustic 



Beyond Caustic 





Two xl 's coalesce to 



Minimiser jumps. 



form a point of inflexion. 
Figure 2: The graph of f( x ,t)(xo) as x crosses the caustic. 



inflexion at = A is the global minimiser. In this case the caustic at x t (X) 
is said to be cool. 

We can also consider the Maxwell set in terms of the reduced action 
function. The Maxwell set corresponds to values of x for which f( x ,t){xo) 
has two critical points at the same action value. If each of this pair of 
critical points also minimises the reduced action, then the inviscid limit of 
the solution to the Burgers equation will jump as shown in Figure [3] and the 
Maxwell set will be described as cool. Note that a Maxwell set can only 
exist in a region with three or more real pre-images if the reduced action is 
continuous. 



Before M t On Cool M t Beyond M t 




Minimiser at x\. Two Xq's at same level. Minimiser jumps. 



Figure 3: The graph of f( x ,t)(^o) as x crosses the Maxwell set. 



We now use the reduced action function to find the Maxwell set. Instead 
of finding the Maxwell set directly it is easier to find the Maxwell-Klein set 
where Definition II .41 is changed to allow x^^xq G C d |22j . 

Theorem 2.8. Let the reduced action function f( x ,t)( x o) be polynomial in all 
space variables. Then the set of possible discontinuities for a d-dimensional 
Burgers fluid velocity field in the inviscid limit is the double discriminant, 

D(t):=D c {D x (/ ( x,t)(A)-c)} = 0, 

where D x (p(x)) is the discriminant of the polynomial p with respect to x. 
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Theorem 2.9. The double discriminant D(t) factorises as, 

D(t) = b 2 m ~ 2 ■ {C t f ■ (B t f 

where B t = is the equation of the Maxwell-Klein set, C t = is the equation 
of the caustic and b is some function oft only; both B t and C t are algebraic 
in x and t. 

Theorem 12.91 gives the Maxwell-Klein set as an algebraic surface. It is 
then necessary to extract the Maxwell set from this by establishing which 
points have real pre-images. Alternatively, we can find the Maxwell set via a 
pre-parameterisation which allows us to restrict the pre-images to be real; for 
this we need the pre-Maxwell set which can also be found using the reduced 
action function. 

Theorem 2.10. The pre-Maxwell set is given by the discriminant, 

{ f(*t(x ),t)( x o) ~ /(*«(*<>),*) O^oA n 

^ { (JTW ) = a 

3 Geometric results 

We now summarise geometric results of DTZ and also [23]. Assume that 
A(x ,x, t) is C 4 in space variables with det (^ g^^ j 7^ 0. 

Definition 3.1. A curve x = £(7), 7 G iV(7o, 5), is said to have a generalised 
cusp at r y = 7o, 7 being an intrinsic variable such as arc length, if £'(70) = 0. 

Lemma 3.2. Let $ t denote the stochastic flow map and $^ 1 r t and Y t be 
some surfaces where if Xq G $^ 1 r t then x = $t(xo) G T t . Then, $ t is a 
differentiate map from &^~ 1 T t to T t with Frechet derivative, 

( d 2 A \ _1 / d 2 A \ 

Lemma 3.3 (2 dims). Let xq(s) be any two dimensional intrinsically parame- 
terised curve, and define x(s) = $i(xo(s)). Let e^ denote the zero eigenvector 
of (d 2 A/(dx ) 2 ) and assume that ker(d 2 A/ (dx ) 2 ) = (e ). Then, there is a 
generalised cusp on x(s) when s = a if and only if either: 

1. there is a generalised cusp on Xq(s) when s = a; or, 

2. Xo(a) G Q^Ct and the tangent -^p(s) at s = a is parallel to eo- 
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Proposition 3.4. The normal to the pre-level surface is, 



n K (x ) 



( d 2 A 



d 2 A 

dxndx 



X(t,x , V5 (i )) 



Proposition 3.5 (2 dims). Assume that atxo G the normal to ^^H^ 

is non-zero, so that the pre-level surface does not have a generalised cusp at 
xq. Then the level surface can only have a generalised cusp at $t(xo) if 
$t(x ) G C t . Moreover, if x = $ t (x ) G ^{^^Ctnf^ 1 ^}, £/ie leve/ 
surface will have a generalised cusp. 

Example 3.6 (The generic Cusp). Figure H] shows how a point lying on three 
level surfaces has three distinct real pre-images each on a separate pre-level 
surface. A cusp only occurs on the corresponding level surface when the 
pre-level surface intersects the pre-caustic. Thus, provided the normal to the 
pre-level surface is well defined, a level surface can only have a cusp on the 
caustic, but it does not have to be cusped when it meets the caustic. 




Figure 4: (a) The pre-level surface (solid line) and pre-caustic (dashed), (b) 
the level surface (solid line) and caustic (dashed), both for the generic Cusp. 



Proposition 3.7. Assume thatx G M t corresponds to exactly two pre-images 
on the pre-Maxwell set, Xq and Xq. Then the normal to the pre-Maxwell set 
at xq is given by, 

f 9 2 A , A / d 2 A . X 1 
n M (x ) = -^— {xo , x ,t)){—{z ,x,t)) ■ 

X(t, Xo, VS (x )) - X(t, Xq, VS (x )) 
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Proposition 3.8 (2 dims). Assume that at x G $7 1 M t; n M (x ) 7^ so that 
the pre-Maxwell set does not have a generalised cusp. Then, the Maxwell set 
can only have a cusp at $t(xo) if <&t{xo) G Ct- Moreover, if x = $t(xo) £ 
<3>t {Q^Ct fl ) ^ e Maxwell set will have a generalised cusp at x. 

Corollary 3.9 (2 dims). Assuming that n-n(xo) 7^ and um{xq) 7^ 0, then 
when the pre-Maxwell set intersects the pre-caustic, it touches a pre-level sur- 
face. Moreover, there is a cusp on the pre-Maxwell set which also intersects 
the same pre-level surface. 

Example 3.10 (The polynomial swallowtail). Let V = 0, k t = and 

So(io,2/o) = x o + x oVo- From Proposition 13. 8\ the cusps on the Maxwell 
set correspond to the intersections of the pre-curves (points 3 and 6 on Fig- 
ure But from Corollary 13. 9\ the cusps on the Maxwell set also correspond 
to the cusps on the pre-Maxwell set (points 2 and 5). Each cusp on the pre- 
Maxwell set lies on the same level surface as a point of intersection between 
the pre-caustic and pre-Maxwell set (ie. 3 and 5, 2 and 6). 




Figure 5: (a) The pre-Maxwell set (solid line) and pre-caustic (dashed), (b) 
the Maxwell set (solid) and caustic (dashed), for the polynomial swallowtail. 

The Maxwell set terminates when it reaches the cusps on the caustic 
(points 1 and 4). These points satisfy the condition for a generalised cusp 
but, instead of appearing cusped, the curve stops and the parameterisation 
begins again in the sense that it maps back exactly onto itself. This follows 
because every point on the Maxwell set has at least two real pre-images, and 
so by pre-parameterising the Maxwell set, we effectively sweep it out twice. 
All of the pre-surfaces touch at the cusps on the caustic. 

These results can be extended to three dimensions. 
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Theorem 3.11 (3 dims). Let 

x e Cus P (H?) = {x e (Q^Ct n fcr 1 ^) = $ t (x ),%(i ) 7^ 0} . 

Then in three dimensions, with probability one, Tn(x) the tangent space to 
the level surface at x is at most one dimensional. 

Theorem 3.12 (3 dims). Let, 

x e Cusp(M t ) = {xe$ t (^Q n $7^) ,x = $ t (x ),n M (x ) ^ 0} . 

Then in three dimensions, with probability one, Tm{x), the tangent space to 
the Maxwell set at x is at most one dimensional. 

4 Swallowtail perestroikas 

The geometry of a caustic or wavefront can suddenly change with singularities 
appearing and disappearing [2J. We consider the formation or collapse of a 
swallowtail using some earlier works of Cayley and Klein. Here we provide a 
summary of results from [22] . 

In Cayley's work on plane algebraic curves, he describes the possible triple 
points of a curve [25] by considering the collapse of systems of double points 
which would lead to the existence of three tangents at a point. The four 
possibilities are shown in Figure El The systems will collapse to form a triple 
point with respectively, three real distinct tangents, three real tangents with 
two coincident, three real tangents all of which are coincident, or one real 
tangent and two complex tangents. We are interested in the interchange 
between the last two cases which Felix Klein investigated [T^ 116] . 




Figure 6: Cayley's triple points. 

As indicated in Section 2, we often parameterise the caustic and level 
surfaces using a pre-parameterisation in which we restrict the parameter to be 
real to only consider points with real pre-images. This does not allow there to 
be any isolated double points on these curves. We now let the pre-parameter 
vary throughout the complex plane and consider when this maps to real 
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points. We begin with a family of curves of the form x t (X) = (x](X),Xt(X)) 
where each x"(A) is real analytic in A. If lm{x t (a + ir/)} = 0, it follows that 
Xt(a + irj) = x t (a — iT]), so this is a " complex double point" of the curve x t (X). 

Proposition 4.1. If a swallowtail on the curve xt(X) collapses to a point 
where A = A when t = t then (A) = -grsKA) = 0. 

Proposition 4.2. Assume that there exists a neighbourhood of A G R sttc/i 
i/iat -^-(A) 7^ fort E (t — 5,t) where S > 0. //a complex double point joins 

the curve x t (X) at A = A u>/ien t = t then -^f(A) = 4xH^) = ^- 

These give a necessary condition for the formation or destruction of a 
swallowtail, and for complex double points to join or leave the main curve. 

Definition 4.3. A family of parameterised curves x t (X), (where A is some 
intrinsic parameter) for which -^f(A) = "rr^PO = is said to have a point 
of swallowtail perestroika when A = A and t = t. 

As with generalised cusps, we have not ruled out further degeneracy at 
these points. Moreover, as Cayley highlighted, these points are not cusped 
and are barely distinguishable from an ordinary point of the curve [25J. 

We now apply these ideas to the caustic where x t (X) will denote the pre- 
parameterisation. The "complex caustic" is found by allowing the parameter 
A to vary over the complex plane. We are interested in the complex double 
points if they join the main caustic at some finite critical time t where rj t — > 
as 1 1 1; at such a point a swallowtail can develop. 

Example 4.4. Let V(x, y) = 0, k t (x, y) = and So(xo, yo) = Xq + xfyo. The 
caustic has no cusps for times t < t and two cusps for times t > t where, 

t = |V2(f)^ = 2.5854... 

At the critical time t the caustic has a point of swallowtail perestroika as 
shown in Figure [71 There are five complex double points before the critical 
time and four afterwards. The remaining complex double points do not join 
the main caustic and so do not influence its behaviour for real times. 




Figure 7: Caustic plotted at corresponding times. 
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Unsurprisingly, these phenomena are not restricted to caustics. There 
is an interplay between the level surfaces and the caustics, characterised by 
their pre-images. 

Proposition 4.5. Assume that in two dimensions at xq 6 ^ 1 H^ H C t 
the normal to the pre-level surface nn(^o) 7^ and the normal to the pre- 
caustic nc{xo) ^ so that the pre-caustic is not cusped at xq. Then nc{xo) 
is parallel to nn(xo) if an d only if there is a generalised cusp on the caustic. 

Corollary 4.6. Assume that in two dimensions at xq G $ < ~ 1 ^ c ^ ®t @t the 
normal to the pre-level surface nu(xo) 7^ 0. Then at <& t (xo) there is a point 
of swallowtail perestroika on the level surface if and only if there is a 
generalised cusp on the caustic C t at <3> t (xo). 

Example 4.7. Let V(x, y) = 0, k t (x,y) = 0, and 5o(xo,2/o) = ^0 + x oVo- 
Consider the behaviour of the level surfaces through a point inside the caustic 
swallowtail at a fixed time as the point is moved through a cusp on the 
caustic. This is illustrated in Figure [HJ Part (a) shows all five of the level 
surfaces through the point demonstrating how three swallowtail level surfaces 
collapse together at the cusp to form a single level surface with a point of 
swallowtail perestroika. Parts (b) and (c) show how one of these swallowtails 
collapses on its own and how its pre-image behaves. 




Figure 8: (a) All level surfaces (solid line) through a point as it crosses the 
caustic (dashed line) at a cusp, (b) one of these level surfaces with its complex 
double point, and (c) its real pre-image. 
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5 Real turbulence 



The geometric results of Section 3 showed that cusps (or in three dimensions 
curves of cusps) on the level surfaces occur where the pre-level surface inter- 
sects the pre-caustic. As time passes, the cusps or curves of cusps will appear 
and disappear on the level surfaces as the pre-curves move. 

Definition 5.1. Real turbulent times are defined to be times t at which there 
exist points where the pre-level surface &~j~ 1 H£ and pre-caustic ^ l C t touch. 

Real turbulent times correspond to times at which there is a change in 
the number of cusps or cusped curves on the level surface iff. 

In (^-dimensions, assuming $ 4 is globally reducible, let /( x ,t)(^o) denote 
the reduced action function and x t (\) denote the pre-parameterisation of the 
caustic. 

Theorem 5.2. The real turbulent times t are given by the zeros of the zeta 
process Q where, 

Ct '■= /(x t (A),t)(Al) - C, 

A satisfies, 

d 

Q^-f(xtW,t)(*i) = fora = l,2,...,d, (5) 
and Xt(X) is on the cool part of the caustic. 

Proof. At real turbulent times there is a change in the cardinality, 

# {A d = Ad(Al, . . . , Arf_i) : /(x- t (Ai,...,A d _ 1 ,A d (Ai,...,A d _ 1 )),t)(Al) = c} . □ 



5.1 White noise in (/-orthogonal directions 



We now consider the Burgers fluid under the potential V{x) =0 and the 
noise Yla=i V k a (x)W a (t) where W a are <i-independent Wiener processes and 
k a {x) = x a with x = (xi, X2, ■ ■ ■ , Xd)- The Burgers equation is then, 



— + (v»-VX = ^Av>*-eW(t), 



(6) 



where W(t) = W 2 (t), W d {t)). 

Proposition 5.3. The stochastic action corresponding to the Burgers equa- 
tion (0|) is, 



A(x ,x,t) 



\x — Xo\ 

2~r~ 



+ |(x 



W(s) ds-ex- W(t) 



\W(s)\ 2 ds + 



2t 



W(s) du 



+ S (x ). 
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Proof. The action is derived as in Section 2 using integration by parts. □ 

Proposition 5.4. If x e t {\) denotes the pre-parameterisation of the random 
caustic for the stochastic Burgers equation (0|) and x° (A) denotes the pre- 
parameterisation of the deterministic caustic (the e = case) then, 

x e t (X) =x°(A) -e / W(u)du. 

Proof. Follows from Definition 11.21 and Theorem I2.2L □ 

Using Propositions 15.31 and 15.41 we can find the zeta process explicitly. 

Theorem 5.5. In d- dimensions, the zeta process for the stochastic Burgers 
equation is, 

Cf = ^o Wjt) (A0-ex?(A)-^ 

where /LJAj) is the deterministic reduced action function, x°(A) is the de- 
terministic caustic and A must satisfy the stochastic equation, 

Vx(f? xUx)j 4\i)-ex t (\)-W(t)) = 0. (7) 

Proof. Follows from Theorem 15.21 having derived by induction the reduced 
action function from Theorem 15.31 □ 

Equation ([7]) shows that the value of A used in the zeta process may be 
either deterministic or random. In the two dimensional case this gives, 

0=(v.yJo w , t) (Ai)-6W(t)).-^(A) > (8) 

which has a deterministic solution for A corresponding to a cusp on the 
deterministic caustic. This point will be returned to in Section 5.3. 

5.2 Recurrence, Strassen and Spitzer 

One of the key properties associated with turbulence is the intermittent re- 
currence of short intervals during which the fluid velocity varies infinitely 
rapidly. 

Using the law of the iterated logarithm, it is a simple matter to show 
formally that if there is a time r such that = 0, then there will be infinitely 
many zeros of Q i n some neighbourhood of r. This will make the set of zeros 
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of Q a perfect set and will result in a short period during which the fluid 
velocity will vary infinitely rapidly. However, this formal argument is not 
rigorous as it will not hold on some set of times t of measure zero [21]. 

The intermittent recurrence of turbulence will be demonstrated if we can 
show that there is an unbounded increasing sequence of times at which the 
zeta process is zero. 

We begin by indicating the derivation of Strassen's form of the law of the 
iterated logarithm from the theory of large deviations [23 [35]. Consider a 
complete separable metric space X with a family of probability measures P e 
defined on the Borel sigma algebra of X. 

Definition 5.6. The family of probability measures P e obeys the large devi- 
ation principle with a rate function I if there exists a lower semicontinuous 
function I : X — > [0, oo] where: 

1. for each I G M the set {x : I(x) < 1} is compact in X , 

2. for each closed set C C X, limsupelnP e (C) < — inf I(x), 

e->o zee 

3. for each open set G C X , liminf elnP e (G) > — inf I(x). 

Let X = Co[0, 1] where C [0, 1] is the space of continuous functions / : 
[0, 1] -> R d with /(0) = 0. Let W(t) be a d-dimensional Wiener process and 
P e be the distribution of ^/eW(t) so that Pi is the Wiener measure. 

Theorem 5.7. For the measure P e the large deviation principle holds with a 
rate function, 

jf n\ _ [ \ Jq f(t) 2 dt : f{t) absolutely continuous and /(0) = 0, 
\ oo : otherwise. 

Definition 5.8. The set of Strassen functions is defined by, 

*r = {/eC [o,i]: 2/(/)<l}. 

Theorem 5.9 (Strassen's Law of the Iterated Logarithm). Let Z n (t) = 
_ i 

(2nlnlnn) 2 W{nt) for n > 2 and < t < 1 where W(t) is a d-dimensional 
Wiener process. For almost all paths u> the subset {Z n (t) : n — 2, 3, . . .} is 
relatively compact with limit set K . 

Following the ideas of RTW, this theorem can be applied to the zeta 
process to demonstrate its recurrence. 
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Corollary 5.10. There exists an unbounded increasing sequence of times t n 
for which Y tn = 0, almost surely, where, 

Y t = W(t)- [ W(s)ds-l [ \W(s)\ 2 ds, 
Jo z Jo 

and W(t) is a d-dimensional Wiener process. 

Proof. If h{n) = (2nlnlnn) _ 5 and x{t) £ K then there exists an increasing 
sequence ni such that, Z ni {t) = h{rii)W{nit) — > x(t), as i — > oo. 

Consider the behaviour of each term in h{ni) 2 n1 l Y t . Firstly, by applying 
Lebesgue's dominated convergence theorem, 

\2 -li 



himyn-'Wim) ■ I W(s)ds -> x(l)- I x(r)dr, (9) 

and, 

himfn; 1 / \W(s)\ 2 ds -> / |x(r)| 2 dr, (10) 



as z — > oo. 



Now let = (x±(t), X2(t), . . . ,Xd(t)) where x a (t) = d~H for each a = 
1,2, ...,d. Therefore, from equations (El) and ( TTOT) . there is an increasing 
sequence of times such that, 



as i — > oo. 

Alternatively, let 



x a (t) 



{ « 2 6 3 



(d)"5* : 0<t<±, 



for a = 1, 2, ... d. Therefore, using equations (jHD and ( fTOl) there is an increas- 
ing sequence of times r, such that, 

v ; 1 1 54 27 18 

Thus, the sequence ti is an unbounded increasing infinite sequence of 
times at which Y t > 0, and the sequence t { is an unbounded increasing 
infinite sequence of times at which Y t < 0. □ 
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d 

Corollary 5.11. If h{t)H~ l ff 0(x) t) {Xi) and h^t' 1 J2 (X) -> 0, f/ien 

' i=0 

iae zeta process Q is recurrent. 

A stronger condition on recurrence can be found in the two dimensional 
case if we work with small e and neglect terms of order e 2 so that, 

C = / ( ,o (A)li )(Ai)-ex°(A)-^(t)- C . (11) 

For this we use Spitzer's theorem, a proof of which can be found in Durrett 

Theorem 5.12 (Spitzer's Theorem). Let D(t) = Di(t) + iD2(t) be a complex 
Brownian motion where D\ and D 2 are independent, -Di(O) = 1 and D 2 (0) = 
0. Define the process 9 t as the continuous process where 6q = and sm(8 t ) = 



pff, . Then, as t — > oo, 

1 [ y dx 



7T 1 + X 2 . 



The process t gives the angle swept out by D(t) in time t, counting 
anti-clockwise loops as — 2tt and clockwise loops as 2tt. 

Let A : M + — > M 2 and consider the behaviour of the process, Y t = A(t) ■ 
W(t). Assuming that A(t) ^ 0, let <pt and 9 t measure the windings around 
the origin of A(t) and W t respectively. Then, Y t = e\A(t)\\W(t)\ cos(4> t - 9 t ). 
Therefore, for Y t = we require cos(0j — 8 t ) = 0, so that the two vectors 
A(t) and W(i) are perpendicular to each other. (Alternatively, this would 
be satisfied trivially if A(t) were periodically zero with t.) 

Corollary 5.13. The small noise zeta process < f77]j is recurrent if there exists 
a bounded function h(t) where h : M + — > M + sitc/i iaai, 



as £ — > oo and iaere exists a function n t such that n t — > oo 

A.n 2 n\ — <\>\ 1 ra t ln£ 

(hit) 2 < 4 an 167r% 2 -40 t 2 + (lnt) 2 ~" °' 

as t — > oo where A{t) = eh(t)x^(\). 
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5.3 Examples in two and three dimensions 

We now consider an explicit example in two dimensions. Since the parameter 
A G K, equation (jSJ) reduces to, 

= ^/(.,(A),)(A)=V,/ ( , t{A) , t) (A)-^(A). (12) 

This gives three different forms of turbulence: 

1. 'zero speed turbulence' where V/^t) (xq) = X(t) = 0. and so the Burg- 
ers fluid has zero velocity. 

2. 'orthogonal turbulence' where V/( It (A),t)(A) is orthogonal to ^f(A) so 
that the caustic tangent is orthogonal to the Burgers fluid velocity. 

3. 'cusped turbulence' where (A) = 0, so there is a generalised cusp on 
the caustic at x t (X). 

As discussed previoiusly, cusped turbulence will occur at deterministic 
values of A and will also correspond to points of swallowtail perestroika on 
the level surfaces. As such, it is not only the simplest form to analyse, 
but also the most important. The categorisation of turbulence leads to a 
factorisation of equation (1121) . 

Example 5.14 (The generic Cusp). For the generic Cusp, the zeta process 
reduces to, 

+e 2 (w(t) ■ jf W{s)ds - ijT \W{s)\ 2 ds 
where A must be a root of, 

= ^\t(2\H 2 - A 2 - 2e{XtW 1 (t) + W 2 (t)}). 

The factor A = corresponds to the cusp on the caustic while the roots of the 
second factor correspond to orthogonal and zero speed turbulence. Firstly, if 
A = then from Corollary 15.111 the zeta process is recurrent. Therefore, the 
turbulence occurring at the cusp on the generic Cusp caustic is recurrent. 

Alternatively, for large times it can be shown formally that the four roots 
which give rise to orthogonal and zero speed turbulence all tend towards 
zero. Thus all four roots tend towards the cusp and consequently, the zeta 
processes associated with each root will be recurrent. 

Moreover, it can be shown that the whole caustic is cool and so all of 
these points of turbulence will be genuine. 
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Next consider the three dimensional case. Thus A e R 2 and equation (JHJ) 
becomes the pair, 

= V,/ (st(A) , t) (A 1 ) ■ ^(A), = VJ {Mm (*i) ■ ^(A). (13) 

In direct correlation to the two dimensional case, we can categorise three 
dimensional turbulence depending on how we solve equations ( I13j) : 

1. 'zero speed turbulence' where again 'V x f( Xt (\),t)(^i) — 0; 

2. 'orthogonal turbulence' where all three vectors ^ x f(x t {X),t)(Xi), af^(A) 
and -ff^(A) are mutually orthogonal. 

3. 'subcaustic turbulence' where the vectors -jf^-(A) and -^f|(A) are linearly 
dependent. The "subcaustic" is the region of the caustic where the 
tangent space drops one or more dimensions. In three dimensions it 
corresponds to folds in the caustic. 

As in the two dimensional case, it follows from Proposition 15.41 that the val- 
ues of A that determine the subcaustic are deterministic. However, unlike 
the two dimensional case, subcaustic turbulence only occurs at points where 
the Burgers fluid velocity is orthogonal to the subcaustic. Hence, we are se- 
lecting random points on a deterministic curve, and so subcaustic turbulence 
involves random values of A. Again the categoristaion of turbulence leads to 
a factorisation in equations ffT3l) . 

Example 5.15 (The butterfly). Let S (x , y , z ) = XqI/o + xlz , this gives a 
butterfly caustic - the three dimensional analogue of the generic Cusp. The 
zeta process is, 

Q = \\\ 2 - ~\\t-4.X\t+-X\Xlt- 12A^A 2 t 2 - 27A^A 2 t 2 + 8A^ 3 
+36 A^ 3 + ^Aft 3 + e((3A 2 A 2 t - AX\t 2 - 9A^)Wi(t) 

-(A 2 + \\t)W 2 (t) + (3AxA 2 + j t - 3X\t - % -X\t)W 3 {t)) 

+e 2 (w(t)- J W(s)ds-^J |H/(s)| 2 ds) -c, 
where A = (Ai, A 2 ) must satisfy, 

= 135t 3 X 9 1 + 96t 3 \l-63\ 2 t 2 \ 6 1 + (16t 3 -8t)\l- (20X 2 + 15eW 1 (t))t 2 \ 4 1 
+ (QX 2 2 - 2 - 6eW 3 (t)) t\l + (A 2 - At 2 eW x {t) - teW 2 (t)) A 2 
+2 (X 2 W x (t) - W 3 {t)) etXy + eX 2 W 3 {t), 

= -27t 2 X 7 1 -12t 2 Xl + 9tXtX 2 + Xl + 3teW 1 (t)X 2 1 + 3eW 3 (t)X 1 -eW 2 (t). 
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Eliminating A2 gives the factorisation, 

= (54t 2 X[ + 6t 2 Xl + Xl + 3teW 1 (t)\l + 3eW 3 (t)X 1 -eW 2 (t)) 

x (Af - 3eW 3 (t)X 1 + 2eW 2 (t)) , (14) 

where the first factor gives zero speed and orthogonal turbulence while the 
second factor gives subcaustic turbulence. 

For large times, it can be shown formally that of the seven roots corre- 
sponding to zero and orthogonal turbulence, five should tend to A = (0, 0) 
and so should give a recurrent zeta process. None of the remaining roots give 
recurrence |21j . 

5.4 The harmonic oscillator potential 

It is not always necessary to resort to Strassen's law to demonstrate the 
recurrence of turbulence; some systems have an inherent periodicity which 
produces such behaviour. The following two dimensional example is taken 
from RTW [M] in which a single Wiener process acts in the x direction. 

Example 5.16. Let k t (x,y) = x, V(x,y) = ^(x 2 u 2 + y 2 uj 2 ) and S (x ,y ) = 
f(x )+g{xo)y where /, /', /'", g, g', g'" are zero when x = a and g"(a) ^ 0. 
Then the zeta process for turbulence at a is given by, 

Q = — - — 2 —^- sin(2u;2t) csc 2 (cj 1 t) {sin{uj\t) f' (a) + U\ cos(cJit)} 2 
4g"{a) 

+ecsc(uJit)Rt — -a 2 Ui sin(2u 1 t) — c, 

where R t is a stochastic process which is well defined for all t. 

Therefore, Q — > ±00 as t — > ^ because csc 2 (/c7r) = 00 where the sign 

depends upon the sign of — & ™i 2 ^^ • Thus, it is possible to construct an 
unbounded increasing sequence of times at which ( t switches between ±00 
and so by continuity and the intermediate value theorem there will almost 
surely exist an increasing unbounded sequence at which Q k =0. 

6 Complex turbulence 

We now consider a completely different approach to turbulence based on 
the work of Section 4. Let (X,Xqq(X)) denote the parameterisation of the 
pre-caustic so that x t (X) = $ t (A,Iq C (A)) is the pre-parameterisation of the 
caustic. When, 

Z t = Im {$<(a + ir],xl c (a + ir/))} , 
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is random, the values of r)(t) for which Z t = will form a stochastic process. 
The zeros of this new process will correspond to points at which the real 
pre-caustic touches the complex pre-caustic. 

Definition 6.1. The complex turbulent times t are defined to be times t when 
the real and complex pre-caustics touch. 

The points at which these surfaces touch correspond to swallowtail pere- 
stroikas on the caustic. 

Theorem 6.2. Let £t(A) denote the pre-parameterisation of the caustic 
where AeM and x t (\) is a real analytic function. If at time t a swallowtail 
perestroika occurs on the caustic when A = A then, 

Assuming that f( x ,t)( x o) is a polynomial in x\ we can use the resultant to 
state explicit conditions for which this holds. 

Lemma 6.3. Let g and h be polynomials of degrees m and n respectively 
with no common roots or zeros. Let f = gh be the product polynomial. Then 
the resultant, 

f m \ n \ fW(Q) \ Ar_1 
where N = m + n and R(g, h) ^ 0. 

Proof. See [22]- □ 

Since f[ Xt r\\ t){ x o) * s a polynomial in xo with real coefficients, its zeros 
are real or occur in complex conjugate pairs. Of the real roots, xq = A is 
repeated. So, 

f(x t (X),t)( X l) = Oo - A) 2 Q(A, i )(4)i/(A,t)(Xo), 

where Q is the product of quadratic factors, 

Q(x,t)(4) = fl{(4-^) 2 + (vl) 2 }, 

8=1 

and Hfx,t)( x o) the product of real factors corresponding to real zeros. This 
gives, 

f(Mx),t)(4)\ xh=x = 2 II {(A - aj) 2 + (vl) 2 } %,*)(A). 

8=1 
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We now assume that the real roots of H are distinct as are the complex 
roots of Q. Denoting fZ t n) t)( x o) ^y ft'W e ^ c ' a simple calculation 

x =X 

gives 

|^A(/f (AU (4) (A))| = 

K t f[(v!) 2 II {K fc - «*) 4 + 2 ((^') 2 + ^) 2 )K fc - <) 2 + (fa?) 2 - W) 2 ) 2 } 

k=l j^k 



x\R x (H,H')\\R x {Q,H)\ 
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i^i being a positive constant. Thus, the condition for a swallowtail pere- 
stroika to occur is that 



P,(t):= i?A(/r(A),A (4) (A)) 



0. 



where we call p v (t) the resultant eta process. 

When the zeros of p v (t) form a perfect set, swallowtails will spontaneously 
appear and disappear on the caustic infinitely rapidly As they do so, the 
geometry of the caustic will rapidly change Moreover, Maxwell sets will be 
created and destroyed with each swallowtail that forms and vanishes as when 
a swallowtail forms it contains a region with two more pre-images than the 
surrounding space. This will add to the turbulent nature of the solution in 
these regions. We call this 'complex turbulence' occurring at the turbulent 
times which are the zeros of the resultant eta process. 

Complex turbulence can be seen as a special case of real turbulence which 
occurs at specific generalised cusps of the caustic. Recall that when a swal- 
lowtail perestroika occurs on a curve, it also satisfies the conditions for having 
a generalised cusp. Thus, the zeros of the resultant eta process must coin- 
cide with some of the zeros of the zeta process for certain forms of cusped 
turbulence. At points where the complex and real pre-caustic touch, the 
real pre-caustic and pre-level surface touch in a particular manner (a double 
touch) since at such a point two swallowtail perestroikas on the level surface 
have coalesced. 

Thus, our separation of complex turbulence from real turbulence can be 
seen as an alternative form of categorisation to that outlined in Section 5.3 
which could be extended to include other perestroikas. 
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